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Mechanical instability takes different forms in various ordered and disordered systems. We study 
the effect of thermal fluctuations in two disordered central-force lattice models near mechanical 
instability: randomly diluted triangular lattice and randomly braced square lattice. These two 
lattices exhibit different scalings for the emergence of rigidity at T = 0 due to their different 
patterns of self stress at the transition. Using analytic theory we show that thermal fluctuations 
stabilize both lattices. In particular, the triangular lattice displays a critical regime in which the 
shear modulus scales as G ~ whereas the square lattice shows G ~ 


PACS numbers: 61.43.-j, 62.20.-x, 46.65.-|-g, 05.70.Jk 

Mechanical instability takes a rich variety of forms in 
systems ranging from granular matter to biological tis¬ 
sue and man-made structures m- Even in lattice mod¬ 
els where there is a remarkably simple principle - the 
Maxwell’s rule which compares discrete numbers of de¬ 
grees of freedom and constraints [T] - multiple classes 
of distinct behaviors arise. Mechanical instabilities in 
lattice models can take the form of either mean-field [la- 
[H] or non mean-field min] transitions, with contin¬ 
uous [n m or discontinuous [T5H22] changes of vari¬ 
ous elastic moduli, and in some case even with a phonon 
structure of nontrivial topology mug. One central con¬ 
cept in understanding these different classes of behaviors 
is self stress, which describes distributions of stress on 
bonds in the lattice with which the net force on every 
site vanishes. This is characterized by the generalized 
Maxwell’s rule [5] 

Nei = dN-Nk + Ns ( 1 ) 

where d is the dimension, N is the number of sites, N}, is 
the number of central force (CF) bonds, Nq is the num¬ 
ber of floppy modes, and Ng is the number of states of 
self stress. According to this rule, each bond added to a 
lattice either removes one floppy mode or increases one 
state of self stress, and in the latter the bond can be con¬ 
sidered as “redundant.” Moreover, states of self stress are 
intimately related to elastic moduli because it describes 
the nature of stress the lattice can sustain [24] . Whether 
a bond is redundant or not strongly depend on the ar¬ 
chitecture of the lattice [53]. Thus a small difference in 
the construction of the lattice can lead to very different 
behaviors near mechanical instability. Different classes of 
lattice models have been found to relate to different ex¬ 
perimental systems, such as jamming of spheres [MIlS], 
glasses 1311], and semi-flexible polymer networks [8]. 

Little is known about how thermal fluctuations af¬ 
fect mechanical instability in these different classes of 
models. In general, fluctuations have been shown to 
change conclusions from mean-field theories about tran¬ 
sitions dramatically, e.g., fluctuations can alter the or¬ 
der and scaling exponents of the transition, or even 
ruling out the ordered phase completely m- It has 


been pointed out that thermal fluctuations may stabi¬ 
lize under-coordinated networks in some special cases 
including crosslinked flexible polymer [SHVldOj . tethered 
networks [STtiSd] . diluted lattices [35l [36], and jammed 
packings [ 571 - 155 ] . but a generic theoretical framework 
for finite-temperature mechanical instability, especially 
in disordered systems, is still lacking. Such generic the¬ 
ory will be important in understanding experimental ob¬ 
servations of mechanical instability which are made at 
finite temperature. 

In this letter, we construct an analytic theory for 
disordered lattices near mechanical instability at finite 
temperature. Our theory treats thermal fluctuations 
through renormalization of elastic moduli and treats lat¬ 
tice disorder through the coherent potential approxima¬ 
tion (CPA) HOI SI]. We use this theory to discuss two 
CF lattice models, the randomly diluted triangular lat¬ 
tice [12] and the randomly braced square lattice [TU] , both 
of which display mechanical instabilities when they pass 
through the CF isostatic point (z) = 2d where z is the 
coordination number. Our main findings are summa¬ 
rized in the phase diagrams in Fig. The triangular 
lattice is stabilized by thermal fluctuations near the iso¬ 
static point, and it displays three regimes of elasticity, 
G ^ (z) — 2d + 0{T) (mechanical), G ~ (critical), 
and G ^ T (entropic), consistent with numerical results 
in Ref. [ 51 ] • Interestingly the square lattice shows very 
different exponents, e.g., G ~ critical regime, 

as well as a nonafhne-affine crossover, originating from 
the sub-extensive number of states of self stress at the 
transition. Our theory is readily generalizable to other 
models to map out finite-temperature phase diagrams for 
different classes of systems near mechanical instability. 

Model — We consider two types of lattices near me¬ 
chanical instability. The first one is the randomly diluted 
CF triangular lattice [Fig. [ija)], in which each nearest- 
neighbor (NN) bond is of spring constant k and is present 
with probability V. Previous numerical studies [13| have 
shown that the T — 0 rigidity transition in the triangular 
lattice is continuous and occur at V — 0.6602 ± 0.0003 
with non-mean-field exponents |43j . The second one is 
the randomly braced CF square lattice [Fig. Bb)] , in 
which all NN bonds are present with spring constant k 


2 



(c) 0.05 
0.04 
0.03 

C~>( 

0.02 

0.01 

0.00 

(d) 0.05 
0.04 


0.55 



G~T 

Entropic 


0.65^ 0.70 

P 








r 0.03 ^ 1 ^! Nonaffine 

^ ^ I Mechanical 

0 / 

/ G-P^k 


0.02 
0.01 


0.2 0.3 

P 


y. 

a 

I 

u 

•p« 

fl 

B ^ 
£ « 


0.4 


0.5 


FIG. 1. Randomly diluted triangular lattice (a), randomly 
braced square lattice (b), with their phase diagrams in the 
space of probability V for having each bond (NNN bonds in 
the case of square lattice) and temperature T measured in 
units of ka^ shown in (c) and (d) respectively. In (d) NNN 
spring constant k = 1. Red dots mark the T = 0 mechanical 
instability point. 


and each next-nearest-neighbor (NNN) bond of spring 
constant k is present with probability V. The T = 0 
rigidity percolation of this model has been studied both 
numerically and analytically, revealing a transition at 
7^ = 0 with mean-field exponents and features of both 
first- and second-order transitions. 

The Hamiltonian of both lattices can be written as 

^ = (2) 

(hi) 


for a deformation that moves positions r) —>■ Ri for each 
lattice site i. Here the sum is over all bonds (*, j), ryij = 1 
if the bond is present and 0 if not, and = G — rj is 
the vector of the bond in the reference space (assuming 
all bonds to be at rest length in the reference space). In 
the triangular lattice kij = k, and in the square lattice 
kij = k for the NN bonds and k^j = k for the NNN 
bonds. 


We choose to use the coherent potential approximation 
(CPA) for the disorder in these models. This method, 
also called the effective medium theory (EMT), has been 
shown to work very well for mechanical properties of dis¬ 
ordered lattices [Mini HOI 11117]. The spirit of the orig¬ 
inal zero temperature CPA is that the disordered lattice 
can be mapped into an effective medium (EM) regular 
lattice in which all bonds are present with an effective 
medium spring constant km as a function of dilution V. 
This mapping is determined by the self-consistency con¬ 
dition that if one randomly chosen bond in the EM is 
replaced according to the rules of the original disordered 
lattice (change into spring constant k with probability V 
and removed with probability 1 — 7^), the average of the 
perturbed phonon Green’s function equals to the unper¬ 
turbed EM Green’s function. 

In this letter we develop a version of the CPA that 
incorporates thermal fluctuation corrections in order to 
study hnite-T mechanical instability. Assuming that the 
mapping to an EM regular lattice still holds at finite tem¬ 
perature, but with a temperature dependent EM spring 
constant fcm(7^,T), we characterize thermal fluctuation 
effects by calculating renormalized elastic moduli of the 
EM as follows. Consider a deformation 

fi ^ Ri = A - fi+ Ui, (3) 

where i labels lattice sites, A denote a macroscopic defor¬ 
mation gradient, and Ui denote finite wavelength fluctu¬ 
ations around A - ri. Eree energy of the lattice at a given 
macroscopic deformation A in finite temperature can be 
obtained by integrating out fluctuations Ui, which yield 
a vibrational entropy term. Although the Hamiltonian 
contains only harmonic springs, expanding H in terms 
of displacements generate higher order terms in displace¬ 
ments, and thus this entropic term depends on A. It 
was shown in Ref. [48] that Ward identity holds for the 
fluctuation corrected free energy F{A). Details of this 
calculation is shown in the Supplementary Information 
(SI). We consider generic deformations of the form 


A = 




( 4 ) 


and taking derivatives of F with respect to t and s yields 
the fluctuation corrected elastic moduli. We then obtain 
renormalized spring constants from these corrected elas¬ 
tic moduli. 

Next we use the CPA procedure to calculate the per¬ 
turbed Green’s function as one bond is replaced 


where in the phonon Green’s function Q renormalized 
spring constants instead of the T — 0 ones are used. Here 
V is the perturbative potential describing changing the 
spring constant of a randomly chosen bond in the EM, 
and 

X = V - V • -V + V ■ • V • - V - ... (6) 
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is the T-matrix of this perturbation. In the case of re¬ 
placing one CF bond this sum for the T-matrix can be 
done exactly. The CPA self-consistency condition can 
then be written in the form 


Tr|fc,=fc + (i-T)r|fc,=o = o. (7) 


We use this equation to determine the EM spring con¬ 
stants. Details are included in the SI. 

Results — We first apply this finite temperature CPA 
to the diluted triangular lattices. The renormalized elas¬ 
tic moduli km is determined by calculating shear rigidity 
of the EM lattice 


I d^F 

1 

(h 

V dt^ 

t=o V8a^/2 

^8 + 2 j 


where a is the lattice constant and ci ~ 0.18. We then 
extract the renormalized spring constant km through the 
relation G = {'/3/4:)km and find 


km — km T 


4ciT 


( 9 ) 


where the last term indicates an increase of EM spring 
constant with T, characterizing fluctuation stabilization 
of this lattice. 

Using this renormalized km we find the CPA self- 
consistency equation 


k 


m 


= 3k 


AT + 


8ciT ■ 
^kmO?- 


( 10 ) 


where AV = V - ^ith = 2/3 being the 

meanfield rigidity transition point for the triangular lat¬ 
tice. Taking T —>■ 0 recovers the zero temperature CPA 
result [33]. In contrast, at finite T, with the addition of 
this fluctuation correction term, we hnd the solution 


k 


m 


3fc|AT| r / 32cir 
2 [ Y 27fca2(AT)2 


( 11 ) 


which has 3 distinct regimes 

[ 3kAV V > and T < 

km = \ ^8cikT/3a^ T > ' , 

ScikT 'p ^ 'pMF,A pnrl 'T 27ka^AV'^ 

\ 9a'^\AV\ f < f c ^ ^ 32ci 

( 12 ) 

as shown in Fig. [^a). Because elastic moduli of the lat¬ 
tice at finite temperature are both proportional to km, 
these 3 regimes correspond to mechanical, critical, and 
entropic regimes of elasticity respectively, as shown in 
Fig. gc). 

We apply a similar calculation to randomly braced 
square lattices. In this case we determine renormalized 
NNN spring constant km by again calculating finite tem¬ 
perature shear modulus, which yields 

(13) 




P 


FIG. 2. (a) CPA solution of km for diluted triangu¬ 
lar lattice as a function of P — Vc- The solutions given 
by Eq. (B14| are shown in solid curves, and the asymp¬ 
totic forms for the mechanical/critical/entropic regimes 
are shown in dotted curve/markers/dashed curves, for 


r/(fcT) = 10” 
CPA solution of 
as a function of V. 




from top to bottom. (b) 
for randomly braced square lattice 
Numerical solutions of Eq. (131 are 
shown in solid curves, and the asymptotic solution for the 
nonafhne mechanical/affine mechanical/critical regimes are 
shown in dotted/dot-dashed/dashed lines, for T/{ka?) = 
10“®, 10“® from top to bottom. 


Comparing to Eq. (B5) we find that the renormaliza¬ 
tion of the NNN bond spring constant contains a non- 
analytic dependence on the ratio kjKm- This comes from 
the unique phonon structure of the square lattice, which 
exhibits lines of states of self stress and floppy phonon 
modes along and qy axis in the first Brillouin zone, 
leading to non-analyticity in phonon Green’s functions. 

The resulting finite-T CPA self-consistency equation 
then reads 


Km 

K 


= v- 



T 

80 ? Km 


(14) 


Its solution can be expanded in 3 regimes too 




( 7 r/ 2 ) 2 fcT 2 X «C 2n‘^ka^V^ and V <C k/tt^ 
Vk T <C 2'K‘^ka^V^ and V 3> k/tt^ 

^ 1/3 2TT^ka^r^ ^T<^ka^ 

(15) 


corresponding to nonafhne mechanical, affine mechani¬ 
cal, and critical regimes, as shown in Eig. [^b) and the 
phase diagram Fig.lud). It is worth noting that the 
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scaling agrees with the ordered lattice results in Ref. [IS] 
with V = Q. 

Discussions — In both lattice models we found that 
thermal fluctuations provide stabilizing effect to elastic¬ 
ity, even when the lattice is below the isostatic point. 
Moreover, the square lattice finite-T phase diagram near 
mechanical instability differs from that of the triangular 
lattice, because of the different structure of states of self 
stress in the square lattice. 

There is an intuitive picture to understand the different 
regimes of the triangular lattice. The mechanical regime 
for V > at low T is straightforward - the lattice 

is stable at T = 0 and thermal fluctuation simple provide 
small correction to elastic moduli. In the critical regime, 
the system starts to have some modes that are of very 
low stiffness. Finite temperature behavior of these modes 
can be understood by considering the following simple 
example [Fig. [^a)]. Here a particle is connected by 2 
springs in the y direction to 2 walls and the springs are at 
rest length when the particle is at a; = 0, ?/ = 0. Assuming 
the 2 springs to be harmonic with force constant k, the 
Hamiltonian of the particle is to leading order 

where £ is the rest length of the springs. It is clear that 
vibrations in x is a floppy mode of this simple system 
because its T = 0 stiffness d'^H/dx^ = 0. Equivalently, 
the “shear modulus” of the system (resistance to trans¬ 
verse displacement t of the top wall) is also 0. At finite 
T, the shear modulus can be calculated using canonical 
ensemble (see SI) 


G = 


d^F 

~W 


C2'/kT 


(17) 



FIG. 3. (a) A schematic picture explaining floppy modes near 
the isostatic point. One particle (red) tethered by two springs, 
showing 0(T^/^) resistance to transverse deformation. One 
example of instantaneous fluctuation configuration is shown 
in gray, (b) A schematic picture of floppy modes further below 
the isostatic point, in which large fluctuations (gray) cost no 
energy. This leads to 0{T) rigidity. 


additional anharmonic terms in the Hamiltonian to pro¬ 
vide stability. An analytic theory of this effect in ordered 
lattices has been discussed in Ref. [45] . 

Our analytic calculations on lattice models provide a 
simple framework to explore fluctuation effects in more 
complicated systems near mechanical instability. For ex¬ 
ample, recent studies on jamming at finite temperature 
revealed significant fluctuation stabilization effect below 
point J {(j) < ft^c) (STHSSj. It would be interesting to con¬ 
sider variations of our lattice models to include effects 
such as one-sided potential and bond-breaking to cap¬ 
ture regimes of elasticity near point J. 


Appendix A: Dynamical matrix and lattice free 
energy 


where C 2 — 0.17. This thermal stiffness for floppy 
modes is an intuitive way to understand the critical 
regime in the diluted triangular lattice, because in this 
regime, there is only a small number of floppy modes Uq 
with positive coefficients for 0{Uq) terms in H, giving 
them entropic rigidity of 0(T^/^). 

When the diluted lattice is deeper below Vc there are 
more floppy modes which allow large fluctuations (flat 
energy landscape, no quartic term), and the correspond¬ 
ing physics, as shown in Fig.3b, is more like a traditional 
“entropic elasticity” picture for flexible polymers. Thus 
the corresponding entropic elasticity is proportional to T 
as in the case of polymers. 

In addition, in these simple pictures and our lattice 
model, at finite T a negative stress is generated by ther¬ 
mal fluctuations, and we take fixed volume ensemble to 
avoid collapsing. However, fluctuation stabilization of 
elasticity is not limited to the case with fixed volume. 
It was shown in Ref. |36j that with zero pressure but 
self-avoidance between bonds, similar stabilization and 
scalings are observed. In our theory, this correspond to 
allow a negative thermal expansion at finite T, and add 


In this section we derive the free energy of an ordered 
lattice at given macroscopic deformations. From the re¬ 
sulting free energy we get fluctuation corrected elastic 
moduli of the effective medium. 

We write the lattice deformation as 

n ^ Rt = A ■ fi + Ui (Al) 

where i labels lattice sites, A is the deformation gradient 
tensor of the (homogeneous) macroscopic deformation, 
and Ui is the fluctuation deviating from the macroscopic 
deformation. To capture bulk and shear moduli we take 

The lattice free energy at given A can be obtained by 
integrating out fluctuations u. To do this we flrst expand 
the Hamiltonian to second order in u 

H = Ho{A) + (A3) 

q 
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where Hq is the energy of the uniformly deformed state, 
Uq is the Fourier transform of fti, and 

D,(A)=u,(A)I + A-M,(A).A^ (A4) 

is the d X d dimensional {d being the spatial dimension) 
dynamical matrix. In this dynamical matrix the scalar 
Vq and second rank tensor are determined by the 
potentials as follows 


Vq{A) 


M,(A) 


= ^ 2[1 - cos{q.lB)] 

B 

= ^2[1-cos(^Zb)] 

B 


Vba 

2\A-IbI 



X 


IbIb 

Ia-ZbP 


Vba \ 
2|A.rBi; 

(A5) 


where is a sum over bonds in one unit cell (here we 
only write down the simple case of one site in each unit 
cell), Ib = Ti — fj is the vector of this bond connecting 
site i and j in the undeformed state, and Vba, ^ba 
represent the bond potential in the uniformly deformed 
state (of deformation gradient A) and its first and second 
derivatives. There is no term linear in Ui in Eq. (A3) be¬ 
cause of the periodicity constraint. A similar derivation 
has been used in Ref. [i5] . 

We then integrate out fluctuations u to obtain lattice 
free energy, 


F{A) = -Tin J 

= ifo(A) + |lndetD,(A). (A6) 

It was shown in Ref. [15] that T(A) only depend on the 
uniform deformation through the combination of strain 
tensor e = (A^ • A — I) so it keeps the same symmetry as 
Hq{A) and thus the Ward identity is satisfied. 


Appendix B: Randomly diluted triangular lattice 


1. The effective medium at finite temperature 


In this section we derive finite temperature elastic 
moduli of the effective medium triangular lattice, and 
obtain fluctuation corrections to the EM spring constant 

■ 

The triangular lattice unit cell contains 3 NN bonds 
with 


h 


a{I,0}, 


h = a 


i 

2’ 2 j 


,h = a 


1 

2’ 2 j’ 
(BI) 


At finite temperature the rigidity of the EM triangular 
lattice in response to shear deformation is given by 


1 d^F 1 2 , ciT\ 

~ VW t=o ~ y/3a^/2 


(B2) 


where the constant 


~ dt^ 


TrlnDAA) 

t=o 


f d^q 52 

J 


trlnD„(A) 

t=o 


~ 0 . 18 , 


(B3) 


where Tr denote a trace of the whole Nd x Nd dimen¬ 
sional matrix (where N is the total number of parti¬ 
cles) and tr denote a trace of the dxd dimensional matrix 
Dq at a given q. This finite temperature shear modulus 
includes contributions from both potential energy in FIq 
and fluctuation correction term which is proportional to 
T. 

Because at T = 0 the shear modulus is solely deter¬ 
mined by spring constant km, at finite T, by identifying 


G-y^k 


(B4) 


we can get from Eq. (B2) the renormalized value of km, 

4ciT 


km — km T 


3a^ 


(B5) 


This equation characterizes fluctuation stabilization of 
the lattice. 


2. The CPA 


In the CPA, the EM spring constant km is determined 
through the following method. We first consider picking 
a random bond between sites j and k from the lattice, 
and replace the spring constant of the bond from km to 
ks- The scattering potential of such a replacement can 
be written as 

= {k, - km){e-^^'^^ - l)(e*« "'= - l)e„e,,, (B6) 

which acts as a perturbation on the dynamical matrix 

D ^ D -k V. (B7) 

The resulting Green’s function can be written in terms 
of a perturbative expansion 

^ = (D-kV)-i (B8) 

where 


where a is the lattice constant. 


g(m) ^ J3-I 


(B9) 
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is the (unperturbed) phonon Green’s function of the EM, 
and 

7" = V - V • • V + V • ■ V • • V - ... 

(BIO) 

is the T-matrix of this perturbation. The CPA self- 
consistency condition that the disorder average of the 
perturbed Green’s function Q to be equal to can 

then be written in the form 

VT\k^=k + {l-r)T\k^=o = 0 . (Bll) 

This equation determines the EM spring constant 
km{V,T). 

In the case of replacing one CE bond this sum for the 
T-matrix can be done exactly, because 


Appendix C: Randomly braced square lattice 

In this model the NN bonds are all present so the NN 
spring constant keeps the same value k, and the CPA 
refers to mapping randomly placed NNN bonds to uni¬ 
form NNN bonds on every site with an EM spring con¬ 
stant Km- 

Each unit cell in the EM has 2 NN bonds and 2 NNN 
bonds, with vectors 

= o{l, 0}, I 2 = a{0,1} 
h = a{l, 1}, h = a{-l, 1}. (Cl) 

In order to extract the renormalized Km we need to 
consider a simple shear of the square lattice, which to 
leading order only deforms the NNN bonds. Similar to 
the calculation for the triangular lattice we have 


V- ■ V- -/ 

'' 1},P '' p,<} 


'■{kg J' 

(2^)12[1 cosip.lB)]eB -Gp ^ • es 

'-{kg 

(B12) 


where es = Ib/\Ib\ is the unit vector along bond T, and 
in the last line we write the integral as h/km where h is 
a dimensionless number. Using this the T-matrix can be 
written as 


T = 



(B13) 


Plugging this into the self-consistency equation (Bll) we 
have 


km _ V -h 
k 1 — h 


(B14) 


At T = 0 we use the bare value km in the EM Green’s 
function and the resulting integral is 

h{0) = 2/3, (B15) 


G = 


1 


V 


t=0 




(C2) 


In contrast to the fluctuation correction in the triangular 
lattice case as in Eq. (B2), here the non-analytic depen¬ 
dence on k and Km comes from the sub-extensive number 
of floppy modes along the Qx and qy directions in mo¬ 
mentum space. These floppy modes provide a divergent 
fluctuation correction to elastic moduli as the lattice ap¬ 
proaches instability Km 0. The integral involving these 
floppy modes has been discussed in detail in Ref. [10]. As 
a result the renormalized NNN spring constant is 


ttT I k 

= + —. (C3) 

Rigorously speaking, the NN spring constant k also get 
renormalized by fluctuation, but because in this problem 
k is always large its renormalization does not have any 
significant effect on stability. 

In the CPA of the randomly braced square lattice we 
replace one arbitrary NNN bond. The calculation for the 
perturbed Green’s function follows a similar procedure 
as we did for the triangular lattice. The self-consistency 
equation takes a similar form 


which give the meanfield T = 0 rigidity percolation point 
for the triangular lattice = 2/3. At T > 0 the 

renormalized value as in Eq. (B5) should be used in the 
Green’s function and we have 


h{T) 


2 kn 

3 ^ 


and the self-consistency equation becomes 




(B16) 


(B17) 


Km _ V - h 
K 1 — h 


(C4) 


with 

h= Km J - cos{p.lB)]eB ■ ■ es (C5) 

where B represent an arbitrary NNN bond. 

The function h for NNN bonds in the square lattice 
behaves rather differently from h in the triangular lattice 
as we discussed above. At T = 0, instead of approaching 
a constant we have 


h = 



which leads to different regimes of entropic elasticity as 
discussed in the main text. 


(G6) 
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with the non-analytic dependence coming from the sub- 
extensive number of floppy modes. This T = 0 solution 
has been discussed in Ref. m- At T > 0, we replace 
Km km in the phonon Green’s function (not in h 
which has an extra prefactor of Km) and expand at low 
temperature. This lead to the finite temperature CPA 
self-consistency equation 



Appendix D: One particle and two springs 


One can then calculate the canonical partition function 
at fixed T by integrating over particle position {x,y), 





dye 




2 



kt'^ \ 

128 € 2 T J ’ 


(D2) 


In this section we discuss the simple picture of one par¬ 
ticle tethered between two walls by two collinear springs. 
To calculate shear modulus, we introduce a transverse 
displacement of the two walls, t/2 on the top and —t/2 
on the bottom. (This is equivalent to fixing one wall and 
only displace the other one by t but that way one need to 
expand around new equilibrium position of the particle.) 
The Hamiltonian is then 


H = 


k 

2 



i^ + t/2)^ \ 
21 ) 



{x — t/2)^ 
2£ 


2 ‘ 


(Dl) 


where K represent modified Bessel function of the second 
kind. Taking derivatives of the resulting free energy F = 
—TlnZ with respect to t gives 


G = 


d^F 

dt"^ 


t=o 


r(3/4) 

2r(i/4) 


0.17V^. 


(D3) 


So the system exhibit a shear modulus that is propor¬ 
tional to \/T. 


[1] J. C. Maxwell, Philos. Mag. 27, 294 (1864). 

[2] C. Calladine, Int. J. Solids Struct. 14, 161 (1978) 

[3] J. C. Phillips, J. Non-Cryst. Solids 34, 153 (1979). 

[4] J. C. Phillips and M. F. Thorpe, Solid State Commun. 
53,699 (1985) 

[5] S. Alexander, Phys. Rep. 296, 65 (1998). 

[6] M. Thorpe, D. Jacobs, M. Chubynsky, and J. Phillips, 

J. Non-Cryst. Solids 266, 859 (2000). 

[7] M. Wyart, Ann. Phys. Fr 30, 1 (2005). 

[8] C. P. Broedersz, X. Mao, T. C. Lubensky, and F. C. 
MacKintosh, Nat. Phys. 7, 983 (2011). 

[9] T. C. Lubensky, C. Kane, X. Mao, A. Souslov, and 

K. Sun, arXiv: 1503.01324 [cond-mat.soft] (2015). 

[10] X. Mao, N. Xu, and T. C. Lubensky, Phys. Rev. Lett. 
104, 085504 (2010) 

[11] W. G. Ellenbroek and X. Mao, Europhys. Lett. 96 (2011). 

[12] X. Mao and T. C. Lubensky, Phys. Rev. E 83, 011111 

( 2011 ) 

[13] D. J. Jacobs and M. F. Thorpe, Phys. Rev. Lett. 75, 4051 
(1995) 

[14] D. J. Jacobs and M. F. Thorpe, Phys. Rev. E 53, 3682 
(1996) 

[15] M. V. Chubynsky and M. F. Thorpe, Phys. Rev. E 76, 
041135 (2007) 

[16] S. P. Obukhov, Phys. Rev. Lett. 74, 4472 (1995) 

[17] C. F. Moukarzel, Phys. Rev. E 68, 056104 (2003). 

[18] S. P. Kasiviswanathan, C. Moore, and L. Theran, in 
SODA, edited by D. Randall (SIAM, 2011) pp. 1237- 
1252. 


[19] O. Rivoire and J. Barre, Phys. Rev. Lett. 97, 148701 
(2006) 

[20] X. Mao, O. Stenull, and T. C. Lubensky, Phys. Rev. E 
87, 042602 (2013) 

[21] J. Barre, M. Lelarge, and D. Mitsche, arXiv:1412.1004vl 
[math.CO] (2014). 

[22] L. Zhang, D. Z. Rocklin, B. G.-g. Chen, and X. Mao, 
arXiv: 1412.6609 [cond-mat.soft] (2014). 

[23] C. L. Kane and T. C. Lubensky, Nat. Phys. 10, 39 (2014). 

[24] S. D. Guest and J. W. Hutchinson, J. Mech. Phys. Solids 
51, 383 (2003). 

[25] K. Sun, A. Souslov, X. Mao, and T. G. Lubensky, Proc. 
Natl. Acad. Sci. U. S. A. 109, 12369 (2012). 

[26] A. J. Liu, S. R. Nagel, W. van Saarloos, and M. Wyart, 
in Dynamical heterogeneities in glasses, colloids, and 
granular media, edited by L. Berthier, G. Biroli, J.-P. 
Bouchaud, L. Gipeletti, and W. van Saarloos (Oxford 
University Press, 2010) Chap. 9. 

[27] J. Cardy, Scaling and renormalization in statistical 
physics (Cambridge University Press, 1996). 

[28] P.-G. De Gennes, Scaling concepts in polymer physics 
(Gornell university press, 1979). 

[29] X. Xing, S. Mukhopadhyay, and P. M. Goldbart, Phys. 
Rev. Lett. 93, 225701 (2004) 

]30] X. Mao, P. M. Goldbart, X. Xing, and A. Zippelius, 
Phys. Rev. E 80, 031140 (2009) 

[31] M. Rubinstein, L. Leibler, and J. Bastide, Phys. Rev. 
Lett. 68, 405 (1992) 

[32] B. Barriere, Journal de Physique I 5, 389 (1995). 













[33] M. J. Bowick, S. M. Catterall, M. Falcioni, G. Thorleifs- 
son, and K. N. Anaenostopoulos, Journal de Physique I 
6, 1321 (1996). 

[34] F. Tessier, D. H. Boal, and D. E. Discher, Phys. Rev. E 
67, 011903 (2003) 

[35] M. Plischke and B. Joos, Phys. Rev. Lett. 80, 4907 
(1998) 

[36] M. Dennison, M. Sheinman, C. Storm, and F. C. MacK- 
intosh, arXiv:1304.3500 [cond-mat.soft] (2013). 

[37] A. Ikeda, L. Berthier, and P. Sollich, Phys. Rev. Lett. 
109, 018301 (2012) 

[38] A. Ikeda, L. Berthier, and G. Biroli, J. Chem. Phys. 138, 
12A507 (2013). 

[39] E. DeGiuli, E. Lerner, and M. Wyart, arXiv:1501.06995 
[cond-mat.soft] (2015). 

[40] P. Soven, Phys. Rev. 178, 1136 (1969) 


[41] E. J. Garboczi and M. F. Thorpe, Phys. Rev. B 31, 7276 
(1985) 

]42] S. Feng and P. N. Sen, Phys. Rev. Lett. 52, 216 (1984) 

[43] These studies were done on generic triangular lattices, in 
which only the topology of the connectivity is kept and 
the periodicity of the lattice is removed to avoid state of 
self stresses associated with parallel bonds. 

[44] S. Feng, M. F. Thorpe, and E. Garboczi, Phys. Rev. B 
31, 276 (1985) 

]45] H. He and M. F. Thorpe, Phys. Rev. Lett. 54, 2107 
(1985) 

]46] W. Tang and M. F. Thorpe, Phys. Rev. B 37, 5539 
(1988) 

]47] X. Mao, O. Stenull, and T. G. Lubensky, Phys. Rev. E 
87, 042601 (2013) 

[48] X. Mao, A. Souslov, G. 1. Mendoza, and T. C. Lubensky, 
Nature Communications 6, 5968 (2015). 






